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Self-similarities and invariant 
densities for model sets 

Michael Baakefl and Robert V. Moodyfl 



ABSTRACT Model sets (also called cut and project sets) are generaliza- 
tions of lattices. Here we show how the self-similarities of model sets are a 
natural replacement for the group of translations of a lattice. This leads us 
to the concept of averaging operators and invariant densities on model sets. 
We prove that invariant densities exist and that they produce absolutely 
continuous invariant measures in internal space. We study the invariant 
densities and their relationships to diffraction, continuous refinement oper- 
ators, and Hutchinson measures. 



1 Model sets and self-similarities 

In this paper we introduce the notion of averaging operators on suitable 
spaces of functions on model sets. An averaging operator encodes informa- 
tion about the entire set of self-similarities with given inflation factor for 
a given model set. It can be interpreted as a Hilbert-Schmidt operator on 
the space of continuous functions on the corresponding acceptance window 
and, remarkably, from this point of view is seen to be an example of the 
recently studied continuous refinement operators. Using this connection we 
can determine the spectrum and associated set of eigenfunctions for any 
inflation factor of any given model set. In particular, the leading eigenvalue 
1 gives rise to an invariant density for the model set. We derive some prop- 
erties of the Bragg spectrum of a model set that has been weighted by an 
invariant density. We also show that an invariant density leads to an ab- 
solutely continuous invariant measure on internal space and we relate this 
measure to a weakly converging sequence of Hutchinson measures. The full 
mathematical development of this work will appear in (2). 
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1.1 Cut and project schemes 

We begin with the notion of a cut and project scheme. By definition, this 
consists of a collection of spaces and mappings: 

R m W n x R™ R" 

U (1.1) 
L 

where R m and R™ are two real spaces, ir 1 and tt 2 are the projection maps 
onto them, and L C K m x R™ is a lattice. We assume that is injective 

and that tt 2 (L) is dense in R™. We call R m (resp. R") the physical (resp. 
internal) space. We will assume that R m and R" are equipped with Eu- 
clidean metrics and that R m x R" is the orthogonal sum of the two spaces. 
For x lying in any of these spaces, \x\ denotes its length. 

A cut and project scheme involves, then, the projection of a lattice into a 
space of smaller dimension, but a lattice that is transversally located with 
respect to the projection maps involved. 

1.2 Example 

A simple, and very useful, example of such a scheme arises from a real 
quadratic irrationality q. We form the ring Z[g] C R and let * be the Z- 
mapping that takes q into its algebraic (quadratic) conjugate. Then the set 
of points Z[q] := {(x,x*) \ x e ^[q]} is a lattice in R 2 and 

R A- R 2 -2^ R 
Z[q] 

where we use the coordinate projections, is a cut a project scheme. An 
important case of this occurs when q = r := (1 + V5 )/2. 

1.3 Model sets 
Let L := tt 1 (L) and let 

()* : L — ► R™ (1.3) 

be the mapping 7T2 ( 7r il£) _1 - This mapping extends naturally to a mapping 
on the rational span QL of L, also denoted by ( )*. Note that the lattice L 
can also be written as 

L = {(x,x*) | x e L}. (1.4) 

Now, let fl C R™. Define 



A = A(O) := {xe L | x* e SI}. 



(1.5) 
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We call such a set A a model set (or cut and project set) if the following 
three conditions are fulfilled, 

Wl ft C 1" is compact. 



W2 ft = hit (ft). 

W3 The boundary of SI has Lebesgue measure 0. 

The mathematical reasons for studying model sets are that they are very 
natural generalizations of lattices, they share many properties with them, 
and they allow symmetries that are otherwise unavailable in lattices of the 
corresponding dimensions. For example, the following properties are shared 
by all model sets A: 

Ml A is uniformly discrete: that is to say, there is an r > so that for all 
distinct x,y G A, \x — y\ > r. 

M2 A is relatively dense: that is to say, there is an R > so that for each 
x € M m the open ball of radius R around x contains a point of A. 

M3 There is a finite set F so that A - A C A - F. 

M4 A - A is a Delone set (i.e. satisfies [Ml, M2]). 

M5 A has a well-defined density d, i.e. 

d = lim — -y^ s . ., = lim — s (1.6) 

s^oo VOl(B s (0)) s^oo Cm s m y ' 

exists, where B s (0) := {x £ W n | |a;| < s} and 

Cm - T (f + 1) (L7) 
is the volume of the unit sphere in W" . 
M6 A diffracts. (See section 3.1 for more on this.) 

A set with the properties Ml and M2 is called a Delone set. A lattice is 
nothing else than a Delone set that is a group. If F = {0} then M3 states 
that A is a g roup, so M3 is in fact a generalization of the group law. The 
limit in ( |l.6|) is easily seen to be independent of the choice of origin for 
the Euclidean space. What is more, it even exists uniformly for sets. This 
means that for any subset S of ft with boundary of measure 0, the relative 
frequency of the points of (A s )* falling into S, as s — > oo, is vol(5)/vol(ft), 
and the convergence is uniform with respect to translation of the set S. For 
more on these properties one may consult || |l3|, [l4|, [llj . 

Model sets arise in situations in which one is looking for Delone structures 
with symmetries that are incompatible with lattices. The most famous 
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example is that of the icosahedral group which cannot appear as the point 
symmetry of any lattice in 3-space. It is known 15, 5| that if G is any finite 
group acting irreducibly in R m and X is any non-trivial orbit of G, then 
either G acts crystallographically in ffi. m , that is to say, there is a G-stable 
lattice of R m ; or there is a G-stable cut and project set in W n that contains 
the set X. This is the origin of the interest in these sets in the theory of 
quasicrystals, see || for background material. 

However, there is a serious price to be paid for moving from lattices to 
model sets. Lattices, by definition, have an entire lattice of translational 
symmet ries. By comparison, in a model set A described by a cut and pr ojec t 
scheme (LI), the set of translational symmetries is the kernel of ()* in (1.3), 
and in all the standard examples this is in fact {0}. Fortunately, in many 
cases of interest, there is nonetheless an abundance of symmetry, as long 
as one is prepared to consider self-similarities instead of group symmetries. 



1.4 Self-similarities 

Definition 1 A self-similarity of A is an ajfine linear mapping t = tq v 

tq v : x h-> Qx + v (1.8) 

on R m that maps A into itself where Q is a (linear) similarity and v € R m . 
Thus Q = qR, i.e. it is made up of an orthogonal transformation R and an 
inflation factor q. 

Let tq v be a self-similarity of A. If A is uniformly discrete, we must 
have \q\ > 1. This is the reason why we also call such a self-similarity an 
affine inflation. We are interested in the entire set of affine inflations with 
the same similarity factor Q. It is convenient to have € A and — 0* £ 
int(f2). Using a translation of A by a vector « gL and the corresponding 
translation of £1 by Vq € L* , we assume this to be the case. This makes no 
structural difference to the set of inflations of A, but simplifies the algebra: 
if G A and tq v is an affine inflation, then v S L and Q(L) C L. In fact, 
we are going to also assume that QL = L. 

Let us then fix once and for all a (linear) similarity transformation Q on 
W 1 such that QL = L. How do we describe the set of all self-similarities 
with similarity factor Q? In preparation for answering this question, it is 
useful to note that there are three different ways of looking at the same cut 
and project set A: first as a Delone set in R m , which we may think of as the 
discrete picture; second as part of the lattice L, which we may think of as 
the arithmetic picture; and finally as a dense subset of Q via the mapping 
()*, which may be thought of as the analytic picture. 

As an illustration of these ideas, note that Q naturally gives rise to an 
automorphism Q of the lattice L, i.e. an element of GL%(L), and a linear 
mapping Q* of M™ that maps O into itself. From the arithmetic nature of 
Q we deduce that the eigenvalues of Q and Q* are algebraic integers and 
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from the compactness of Q that Q* is contractive. Furthermore, one can 
deduce || that Q* is diagonalizable from the corresponding property of 
Q. Strictly speaking we can only deduce that the eigenvalues of Q* do not 
exceed 1 in absolute value, but we will always assume that in fact they are 
less than one in absolute value. In the sequel we will normally denote the 
contraction Q* by A to match various sources we will refer to frequently. 
Define 

tt Q := {jiel"|ifi + «cO}, (1.9) 

we say that Q is compatible with A if int(f2g) =^ 0. In this article, we shall 
always assume that not only f2, but also fig is Riemann measurable, i.e. 
dflg has zero Lebesgue measure. Interpreting ( JL.8[ ) on the window side we 
obtain: 

Proposition 1 Let A = A(f2) be a model set based on a window Q that 
satisfies the window conditions Wl - W3. Let Q be a similarity compatible 
with A. Then the set Tq of affine inflations with the same similarity Q is 
the set of mappings tq v where v runs through the set 

T Q = {veL\v* en Q }. (l.io) 

In particular, Tq is also a model set. 

Let us pause to consider the special situation where Q is convex. In this 
case, rig is also convex and hence satisfies the conditions that we need. 
If in addition Q* — e ■ t, < e < 1, which actually is often the case in 
examples of physical relevance, one obtains 

n Q = (l-e)fi. (1.11) 



If -1 < e < 0, but fl = -O, (|l.llp is still true if e is replaced by |e|. This 
happens in our Examples. 

1.5 Example 

In Example [b^ above, take q = r. A simple model set is defined by 

A = {xe Z[t] I x e [-1,1]}. (1.12) 

Let us look at the inflation factor q — r. Multiplication by q determines 
the contraction A = r'-l, r' = — 1/r, on the internal side and 

n T = {u e R I -- ■ [-1, lliuc [-1, i]} = [-1 + - , l - -] . (i.i3) 

r r r 

Thus, for all v E Z[r] with v' E [-1/r 2 , 1/r 2 ], 

t TV : x i— > ra; + w (1-14) 
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is an inflation of A, and 

*;,„ : V» —+v' (1-15) 

T 

is the corresponding contraction in internal space. 

2 Averaging operators and invariant densities 

One of the most famliar techniques in the theory of group representations 
is the use of group averages in order to produce invariants. Thus for a finite 
group G one typically invokes the process 

geG 

which averages the function F over the group, where g ■ F(x) := F(g x). 
We intend to do exactly the same thing replacing G by the set of all self- 
similarities Tq of a model set A. Since Tq will be infinite we have to be a 
little careful in averaging. 

For any subset T C M m and for any s > 0, we thus define 

T s := {x e T | |x| < s} (1.17) 

where \x\ = (x ■ x) 1 / 2 is the standard Euclidean norm on R m . 



2.1 Averaging operators 

Definition 2 Let A = A(f2) be a model set based on a window Q that 
satisfies the window conditions Wl - W3. Let Q be a similarity (with 
inflation factor q) that is compatible with A, let T = Tq be the set of all 
self- similarities of A with similarity Q, and let T = Tq be the corresponding 
set of translations. Let p : L — > R be any function on L that vanishes off 
A. Then the average of p over T is defined to be 

provided this limit exists. We say that p is a Q-invariant density on A if 

ID1 p is non-negative on A, 

ID2 Ap = p, and 

ID3 p is normalized, i.e. 



(1.19) 
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Note that for x <E A and v € T, t v l x € L but t v x does not in general 



lie in A. Thus in (1.18) one can expect that many of the summands on the 
right hand side will be because p vanishes off A. Note also that ( L.lSj ) is a 



normalization per point of A. This can be changed to a normalization per 
unit volume if necessary because the density d of points of A exists. 
Let C(fi) be the space of all continuous complex- valued functions on 



with support in Q. Via the mapping ()* of (1.3) we obtain a space C(A) of 
functions on L, vanishing off A: for / € C(f2), we define p = p^ E C(A) by 

p(x) = vol(n) •/(**) , (1.20) 

for x £ A, and p{x) = otherwise. Here, the normalization constant vol(fi) 
is thrown in to make things more convenient later on. 



Bearing in mind (1.20), we rewrite (1.18) as a new averaging operator 
on C{9): 

|det(Q)| 



(Af)(x*) = lim J2 /((O-V) (1.21) 



v£T s 

for all x* £ A*. 

Now, it is well-known (l^, |l7|, || that the point s 7r 2 (L) are uniformly 
distributed in R™, in the sense described in Section |l.3| . In particular, the 
points of T* are uniformly distributed in f2g, hence lim^oo (T s )* = Qq. 
Using Weyl's theorem [[l8| [l2|, the continuity of /, and the fact that £Iq is 
Riemann integrable, we obtain 

Using the fact that det(Q) = det(A) -1 (remember that A = Q*) and 
introducing the normalized indicator (or characteristic) functions 

Xs := (1-23) 
vol(b) 



defined for all measurable subsets S of R™, we can rewrite ( 1.22 ) in a 
number of equivalent ways: 



(Af)(x) = I X Ay+nQ {x)f{y)dy 
1 



\&et{A)\ 
1 



X A -i nQ (A i x-y)f(y)dy 

|det(A)| , X nQ (x-y)f(A^y)dy. (1.24) 

This shows that the averaging operator A is a continuous refinement op 
erator in the sense of || in the one-dimensional case and |l(| in the multi- 
dimensional case. In Fourier space, by application of the convolution theo 



rem, (1.24) reads as 

(Af)(k) = X^WftA'k) (1.25) 
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where f(k) := J R „ e~ 2 * lk - x f{x)dx and f(x) = L n e 27 "^ f(k)dk. 



2.2 Invariant densities 

At this point, we can determine a function f — f p corresponding to an 
invariant density p. We are looking for a 1-eigenfunction of the operator 
A. The normalization condition ID3, seen on the w indow side, becomes, 
using Weyl again, L„ f(u)du = 1. Iteration of ( 1.25| ) leads to 



/» = Xn Q (k)-Xn Q (A t k).....X nQ ((A t ) N k)-f((A t ) N+1 k). (1.26) 

Since A is a contraction, so is A*. Consequently, (A t ) N+1 k — > as iV — > oo. 
Since / and -X"n Q are C°° and /(0) = 1, we can take the limit and obtain 

/» = n^ Q ((^) = n : 0l( o ) ■ w 

N=0 N=0 y Q> 

This is an infinite product with compact convergence that solves our prob- 
lem: / is an infinite product of C°°-functions, and is itself C°°. The func- 
tion / is now the inverse Fourier transform of /. By construction, it is the 
Radon-Nikodym derivative, and hence the density, of an absolutely con- 



tinuous invariant ( 1.43 ) measure (with respect to Lebesgue measure) in 



internal space. Again applying the convolution theorem we arrive at: 

Proposition 2 Let A = A(f2) be a model set based on a window f2 that 
satisfies the window conditions Wl - W3. Let Q be a similarity compatible 
with A and let A := Q* . Then there is a unique Q-invariant density p for 



A lying in C(A). This is given through p — pp see (1.20) 7 where f is given 
by the infinite convolution product 

OO 1 

^^voip^y- (1 - 28) 

Note that this convolution of characteristic functions defines a C°° function 
with compact support contained in fi. 

If Q,Q 2 ,..., are all compatible with A then it is instructive to look 
at the corresponding invariant densities / = fn\, /(2), ■ ■ ■ The sets {fig™} 
are increasing and (J VLq™ = Q. The functions of the sequence {/(&)} 
become increasingly concentrated around and it is natural to expect 
linis^oo // s ) = 1q/vo1(£1). This is illustrated in the case of our Example 
by the sequence of graphs of Fig.l. 



2.3 Example 



We continue with Example 1.5 from above and determine the invariant 
density (or rather its Fourier transform) on [— 1 , 1] corresponding to the 
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inflation factor t. If / is this density, then equation (1.25) becomes 



/(*) = yi[-i/TM/r=](*0/(-W 

Routine calculation gives: 

sin(27rafc) 



dx = 



So, we obtain 



f(k) 



sin 



2 ^/(- fc /-) 



(^) 



n / 2 / 

AT=2 I T7f 



(1.29) 



(1.30) 



(1.31) 



The calculations for other inflation factors are similar; Fig. 1 shows some 
of the resulting invariant densities. 
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FIGURE 1. Invariant densities for the model set of Example 
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for inflation 



2.4 Eigenvalues and eigenf unctions 

According to JTo| , if the eigenvalues of A* are {a 1 , . . . , a n }, < |a,-| < 1, 
then the spectrum of A is 



spec(^) = {a a I a £ Z^ } 



(1.32) 
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where a a := a" 1 • . . . • a® n . Furthermore, the multiplicities are those sug- 
gested by the notation: 

mult(A) = #{a | a a = A} . (1.33) 
In particular, if a x = . . . = a„ = a, then A = a' a ' , \a\ = a 1 + . . . + a n , and 



mult (A) = 



\a\ +n- 1 
lal 



a formula that is well-known from the rt-dimensional harmonic oscillator. 

Remark: In |], [l(]] A is assumed to be diagonalizable and all the eigenval- 
ues are assumed real. The diagonalizability of Q guarantees that of Q* = A 
PJ. By allowing complex- valued functions to enter the picture, it is not 
necessary to assume that the eigenvalues of A are real. Nor do we wish to 
impose such an assumption since we do not have this type of control over 
the spectrum of A. 

It is quite easy to find eigenfunctions representing these eigenvalues 
in terms of the invariant density /. To do so, we choose an eigenbasis 
{vi, . . . ,v n } for A* in C™, using the fact that A* is diagonalizable. Any 
k € C" can be written as k — k 1 v 1 + . . . + n n v n and the Kj = Kj(k) are the 
corresponding coordinate functions (which we allow to be C- valued). Now, 
fix b e Z™ , define n b = k\ x ■ . . . ■ and let u(k) := (n(k)) b f (k) . Then 



we obtain from ( 1.25 ) 



Au(k) = Xa Q (k)u{A*k) - X^k^niA^ffiA'k) 

= a b (n(k)) b f(k) = a b u(k). (1.34) 

Returning from the Fourier domain, this enables us to write down the 
eigenfunctions of A. If {v\, . . . , v^} denotes the dual basis (i.e. vj-Vj — <J..), 

we define the directional derivative T)j := uj • V and obtain 
Proposition 3 The partial derivatives 

D b f, b e Zl Q (1.35) 

are eigenfunctions of the refinement operator A, with eigenvalue a b . 

Some of these derivatives, calculated for our guiding Example, are shown 
in Figure 2. 



3 Further remarks 
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FIGURE 2. The first 4 derivatives of the invariant density of Example 
inflation factor t. 

3. 1 Diffraction and a product formula 



1.5 



for 



The product formula ( 1.27j ) has an interpretation on the physical side of 
our picture. To see this, we define functions g s and h s by 



— E 



27zik-v 



9s{k) 



(1.36) 



-2irik-' 



weA„ 



and, furthermore, functions g and h by 

g{k) = lim g s (k) 

s — ► OC 

h(k) = lim h s {k) 



(1.37) 



The existence of g as a well-defined function on R m is a known conse- 
quence of the fact that T is a model set. We now have [Q: 

oc 

Proposition 4 fe(fc) = fl odQ^k). 

N=0 

The significance of the functions g and /i appears in the context of diffrac- 
tion. Suppose that w : L — ► R>o is some bounded non-negative function 
that vanishes off A. Define the tempered distribution 



(1.38) 



xeA 
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where S x is the Dirac measure at x. The limit, as s — > oo, of the averaged 
auto-correlation of this measure, which exists for model sets, is the auto- 
correlation measure of A (also called its Patterson function, though it is a 
distribution) 



Its Fourier transform is a positive measure -y w which is the diffraction pat- 
tern of A. The point part of this measure is the Bragg spectrum of A. In the 
case that w is the indicator function (i.e constant value 1) on a model set 
A then ^ w is a pure point measure and we say that A has a pure point spec- 
trum. In any case, the Bragg spectrum can be calculated from the simpler 
function 



Our functions g and h correspond to the cases when w — 1\ and w is 
the invariant density p, respectively. In particular, h allows us an exact 
description of the intensities of the Bragg spectrum of (A,p). Notice from 
the product formula that its support necessarily lies inside the support of 
the Bragg spectrum of A itself. It goes without saying that both g and h 
are highly discontinuous functions which, for a model set, are non-zero only 
on a dense point set of zero Lebesgue measure. This set is contained in the 
so-called Fourier module of the model set, i.e. the in the set ir 1 (L°), where 
L° = {ye R' m +™ I x ■ y e Z for all x € L} is the dual lattice of L. 

3.2 Hutchinson measures 

Let / be the invariant density of C(f2) corresponding to the compatible 
similarity Q on A = A(Q). There is a corresponding measure /u = \ip with 
support contained in Q, defined by 




(1.39) 




(1.40) 



provided that this limit exits everywhere. In fact [|6| 

\g(k;w)\ 2 = ■%,({*}), for all fee If 



(1.41) 




(1.42) 



We have n f (Sl) = 1. The : 
define t* v ■ fi f by t* ■ fJL f {Y) 





(1.43) 
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Now fix some s > and consider the finite set of contractions t* on f2, 
where v € T s . According to there is a unique non-negative Borel measure 
/i s on Q for which /U s (f2) = 1 and which is invariant in the sense that 



1 

" vGT 3 



Furthermore, this measure is the unique fixed point of the process on the 
set of regular Borel measures on Q with mass 1 that averages a measure 
by the right hand side of (1.44). We call this a Hutchinson measure. Using 
the Levy continuity theorem one can show 

Proposition 5 The sequence of Hutchinson measures fulfils: 

(i) fi s (k) = UN=o9.((A t ) N k), for any s > 0; 

(ii) {fj. s } ~ * fi-fy where the convergence is uniform on compact sets; 

(iii) {fi s } — ► M/? i n the sense of weak convergence, i.e. {fi s (ip)} — > Hf(<p) 
for all f € C(f2). 

In particular, for any function ip S C(f2), fi s {ip), for large s, is a good 
estimate of jxAip). Also, starting from any probability measure on f2, the 
Hutchinson iteration of (1.44) will converge to fi 8 and this procedure, for 
large s, will also give a good approximation to /i^. 



3.3 The topology of L and A 

The space of functions C(A) on which the averaging operator ^4 acts seems 
strange. It is defined in terms of the topology of Q, in R™ and this is very 
different from the discrete topology that we see on A induced by the topol- 
ogy of its ambient space M. m . The appropriate topology for A (and L) is 
defined intrinsically as follows O] . For each compact set K of K. m define 

N K = N K (A) := {v e L \ v + (An K) = An (v + K)}. (1.45) 

Thus Nk is the set of vectors v for which translation by v is a bijection of 
the if-patch of A and onto the (v + -ftT)-patch of A. Note that the mapping 
K i— * Nk is inclusion reversing. 

Proposition 6 |l(| |l7| Suppose that L* n dQ = 0. Then the collection 
of sets {Nk \ K C M m , K compact} is a basis of neighbourhoods of 
for a topological group structure on L. Furthermore, ffi™ with its standard 
topology is the completion of L under the mapping ()* : L — > W 1 . With this 
topology on L, the space C(A) is precisely the space of continuous functions 
on L whose support lies in A. 

The intuition behind continuity of a function <f> (defined on L) with respect 
to this topology is this: if translation by v is a bijection of two "large" 
patches of A, then (f>(x + v) — 4>{x) is (uniformly) small. 
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4 Outlook 



The existence of positive invariant densities naturally suggests probabilistic 
interpretations. In || we study the spectral properties of certain stochastic 
sets whose sites are selected from that of a model set on a probabilistic 
basis according to the density p. Effectively, this gives a distribution of 
points which, after *-mapping them into internal space, looks like our win- 
dow shaped by the invariant density /. This may provide an alternative 
explanation of the recently made observations of such profiles in real data 
fr) . Furthermore, as such sets do have finite entropy density, they might 
be useful for further models of cntropic stabilization of quasicrystals. 

In this article, we have focused on one-component model sets. It is im- 
portant for multi-component or multi-coloured systems to be able to realize 
the similarity-averaging process and the existence of invariant densities in a 
matrix generalization of what we have done here. This means that we have 
a finite family of model sets based on cosets of a common Z-module and 
matrices of similarity maps between them. In fact this set-up results in ma- 
trix continuous refinement operators and ultimately again in the existence 
of invariant densities. What is particularly interesting in this case is the 
appearance of a Markov matrix of weights for the contributions from the 
various windows relative to each other. An exposition of this will appear 
in §. 
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